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ABSTRACT 

We implement the efBcient line of sight method to calculate the anisotropy and polarization of the 
cosmic microwave background for scalar and tensor modes in almost- Friedmann- Robertson- Walker mod- 
els with positive spatial curvature. We present new results for the polarization power spectra in such 
models. 

Subject headings: cosmic microwave background — cosmology: theory 



L INTRODUCTION 



— The aniootropy of the coomic microwave background 
(CMB) plays a key role in many areas of modern cosmol- 
ogy. Joint anal yses of current CMB and type la s uper- 



novae data (e.g. Efstathiou ct al. 199E | ; T urner 1999 ^ 



goot that the univeroo io within a factor fitot — I 0^ = 
1 ± 0.2 of the critical density required for a flat geometry. 
Closed models (in which fitot > 1) therefore account for an 
important sector of the possible parameter space. More- 
over, maximum likelihood searches require theoretical pre- 
dictions over a much larger volume of parameter space to 
establish reliable error estimates on the parameters under 
consideration. It is therefore vital to have a fast and ac- 
curate method for calculating anisotropies for models at 
least within the range 0.4 < 17tnt < 1-6. Previous pa rame- 
ter determinations, such as lEfstathiou et al. (1999 ), have 



had to rely on analytical approximations for rapid calcu- 
lations of the CMB power spectrum in the closed region, 
since the current st ate of the art codes, such as the widely- 
used CMBFAST ( ^eljak fc Zaldarriaga 199^ , do not yet 
support closed models. 

In this Letter we describe a numerical implementation 
of the linearized equat ions of the 1+3 covariant approach 



to CMB a. nis otropies (Challinor fc Lasenby 1999 ; Challi- 
nor 19994 H; pebbie, Dunsby, fc Ellis 1999| ) in almost- 



Friedmann-Robertson- Walker (FRW) models with open, 
flat, and closed background geometries. Our code thus 
allows an efficient exploration of the full cosmological pa- 
rameter space. We present new results for the polariza- 
tion power spectra from scalar and tensor perturbations 
in closed models. The intensity po wer spectra in closed 
models have been calculated before: White fc Scott (1996 ) 
integrated the full Boltzmann hierarchy directly for scalar 
perturb ations, extending the earli er semi -analytic predic- 
tions of lAbbott fc Schaefer (19861 ). Also, [Allen, Caldwellj 
fc Koranda (1995| ) used the semi-analytic approach (ad- 



equate on large scales) to calculate the tensor spectrum. 
None of these calculations included the effects of polariza- 
tion. 

The 1-1-3 covariant formalism provides a physically 
transparent, exact (fully non- linear) description of both 
dynamics and radiative transfer in general cosmological 
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mode ls ( [Ellis fc van Elst 1999| ; [Maartens, Gebbie, fc Elhs 
1999). The full formalism admits a gauge- invariant lin- 



earization about FRW models, resulting in a linear per- 
turbation theory which is arguably simpler, and more 
physically transpa. r ent. than other approaches (e.g. Ma fc 



Bertschinger 1995[ ; [Hu et al. 1998| ). Our implementation 



of the 1+3 formalism is based on the field-tested CMB- 



FAST code written by [Seljak fc Zaldarriaga (1996[ ). Their 
code uses a line of sight integration method to achieve high 
efficiency without compromising accuracy. 

2. IMPLEMENTATION 

2.1. Basic formalism 

We employ the 1-1-3 covariant approach to perturbations 
in cosmology (e.g. Ellis fc van Elst 1999 ), in which depar- 
tures from exact FRW symmetry are described by gauge- 
invariant variables derived from physical observables rela- 
tive to some timelike 4- velocity field The equations of 
radiative transfer can be recast as propagation equations 
along the integral curves of u'^ for the multipoles of the 
total intensity, Ia^ , and the electric and magnetic compo- 
nents of th e total linear polarization, and Ba^ ( [Chalfi- 



|nor 1999a[ ). Here the multipoles are projected (relative to 
u"-) symmetric trace free (PSTF) tensors {Ai represents 
the index string oi . . . ai), which provide a basis-free alter- 
native to the more common scalar-valued multipole coef- 
ficients in spheric al harmonic expansions of the intensity 
and polarization (Kamionkowski, Kosowsky, fc Stebbing 
1997; Seljak fc Zaldarriaga 1997 ). We only consider linear 
polarization here, since circular polarization is not gener- 
ated by Thomson scattering. For small departures from an 
FRW model, the intensity multipoles evolve as (Thornc 



1981 



Lasenby 199£) 



EUis, Matravers, fc Treciokas 198^ ; [Challinor & 
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CMB anisotropies in closed FRW models 



where Q is the expansion of u", a at is the shear, and Aa is 
the acceleration. Here, D" is the totally projected covari- 
ant derivative, an overdot denotes the covariant derivative 
along and angle brackets denote the PSTF part of the 



enclosed indices. The electron number density is no in its 
rest frame which has relative velocity w", and the Thom- 
son c ross section is (Tt - For the electric polarization, we 
have dChallinor 1999a|) 



' 3 ' (i + l)^ — 



-ncCTT 
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and for the magnetic polarization: 



bA, 



£ 2 

Dt^at^Ae-i) i , -i^CUrlgAt = ^icP-T^At ■ (3) 



(2£+l) 



(^ + 1) 



The evolution of the electric and magnetic multipoles 
are coupled through the curl terms, where curlf^^ = 
€bc{ai'D^£Ai^i)'^\ with Cabc the alternating tensor in the rest 
space of u"'. The multipole equations hold for a gen- 

eral linear perturbation around an arbitrary FRW model, 
and for a general choice of u". The equations must be 
supplemented by the cova riant hydrodynamic and gravi- 
todynamic equations (e.g. Ellis fc van Elst 1999 ) to de- 
termine the first-order source terms (Tq6, Aa and Va- For 
an alternative approach to polari zed radiative tr ansfer in 
general FRW geometries, see e.g. |Hu et al. (1991 ). 

The dimensionless power spectrum of the i ntensity 
aniso tropies i s defined by the ensemble average ( iGebbiel 
& EUjs 1998|; [Challinor fc Lasenby 19991) 



in 



{2£y. 



{Air i-meiy 



:{IaJ' 



(4) 



Power spectra for the electric polarization multipoles, ^, 
the magnetic multipoles, Cf^, and the cross-correlation 
between electric polariza tion and the inte nsity, C^^, can 
be d efined analogously (Ghallinor 1999a ). (To conform 
with Seljak fc Zaldarriaga 1997| we include a factor of 
v/[(£ + 1)(£ -I- 2)/l[i - 1)J on the right-hand side of eq. [| 
for each factor of the polarization.) 

We solve the multipole equations (0-||) by expanding 
the first-order variables in PSTF tensors derived from the 
appropriate scalar, vector, and tensor eigenfunctions of the 
comoving Laplacian 5^D°Dq, where S is the scale factor, 
with eigenvalue . The different perturbation types de- 
couple at linear order, with each giving rise to a set of 
coupled first-order differential equations. The Boltzmann 
multipole e quations for scalar and tensor modes are given 
in detail in Challinor (1999c). The equations describing 
perturbations in the other matter components and the ge- 
ometry, which determine the sourc e terms in the Boltz- 



mann hierarchies, can be found in Challinor fc Lasenby 



(1999§; pm"llinor (1999h| ); [Gebbie et al. (19991 ). The mode 
expanded Boltzmann equations can be solved formally as 



integrals along the line of sight. These integral solutions 
form the basis of the line of sight algorithm employed 
by CMB FAST. The solutions for t he intensity are given 
in the 1-1-3 covariant formalism in Challinor fc Lasenby 
(199S| ); [Challinor (1999b| ); [Gcbbic et al. (1999|); solutions 



for the polarization are given in Challinor (1999c). Equiv 



alent res ults in the total angular momentum formalism are 
given in Hu et al. (199S). The present day multipoles are 



obtained by integrating the product of source functions, 
which are inexpensive to compute, with special functions 
(derived from the hyperspherical Bessel functions) which 
result from the projections of the eigenfunctions of the 
Laplacian onto directions on the sky. 

For scalar modes, define = (fc^ -I- K)/\K\, where 
&K/S'^ is the curvature scalar of the spatial sections in 
the background FRW model. In the closed case, the reg- 
ular scalar eigenfunctions of the comoving Laplacian are 



Abbott fc 



complete for v an integer > 1 (e.g. [Tomita 198: 

Schaefer 19861 ). The mode with v — 1 does not contribute 



^ap 

i 



to the perturbations, while the modes with v — 2 (which 



can o nly represent isocurvature perturbations; Bardeen 



1980) only contribute to the CMB dipole. For tensor 
modes = (JT + 'iK) /\K\, and the regular eigentensors of 
the Laplacian are complete with v an integer > 3; ex plicit 
PSTF representations are given in Challinor (19991:). 



2.2. Numerical Evaluation 



The implementati on strategy of CMBFAST ( [Seljak 
fc Zaldarriaga 1996 ) requires only minor modifications 
for closed models. The modes are now discretized with 
wavenumber > 3. A given v only contributes to multi- 
poles with £ < v, so, unlike the open case, the Boltzmann 
hierarchies for a given i' truncate at finite £. For large ly, 
it is possible to terminate the hierarchies artificially (tak- 
ing care to avoid spurious reflection of power) at some 
lower £ without compromising accuracy in the evaluation 
of the source functions. The sources are calculated at ap- 
proximately logarithmically spaced integer values of v and 
interpolated as needed. 

In closed models, a given linear scale at last scattering 
subtends a larger angle on the sky today than in open 
or flat models. This geometric effect shifts power in the 
CMB spectra to smaller £. Since CMBFAST only com- 
putes the power spectra at a few values of £, we adjust 
the ^-sampling according to the the curvature to maintain 
accurate interpolation in all cases. 

For the line of sight integral over sources, we calculate 
the hyperspherical Bessel functions by integrating a sec- 
ond o rder differential equation, as in the open C MBFAST 
code ( Zaldarriaga, Seljak, fc Bertschinger 199g ) The start- 



ing values for the Bessel function and its derivative are 
found using a recursive evaluation or, for the cas es where it 
is acc urate and faster, the WKB approximation (Kosowsky 
1998 ). We assume that the development angle x satisfies 
X < TT, where x — \/\K\rio and 770 is the conformal age of 
the universe. This permits calculations with models up to 
f^tot ~ 1-7 for a matter fraction flm ~ 0.4. For the case 
X ~ TT, all lines of sight converge to the antipodal point 
close to last scattering, and the power spectrum becomes 
featureless ( [White fc Scott 1996[ ). 

It is important to maintain numerical stability in the 
differential equations for the dependent variables. This 
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is especially true for scalar perturbations with isocurva- 
ture initial conditions, where a poor choice of dependent 
variables can lead to large violations of the Einstein con- 
straint equations unless the initial conditions are specified 
to exquisite accuracy. (This problem is particularly acute 
in the Newtonian gauge, e.g. |Ma fc Bertschinger 1995| .) We 
choose to work in a frame in which coincides with the 4- 
velocity of the cold dark matter (CDM). For scalar pertur- 
bations, we determine the perturbations to the geometry 
by evolving the projected gradient of the 3-curvature on 
hypersurfaces orthogonal to which gives good numer- 
ical stability for all initial conditions. (With this choice 
of dependent variables, our equation set is equivalent to a 
subset of the widely-used synchronous equations, gauge- 
fixed to the CDM.) We allow for adiabatic and two types 
of isocurvature initial conditions for scalar perturbations. 
Also, we provide the option of specifying the primordial 
power spectra in non-parametric form, which is useful for 
inverse problems such as initial power spectrum recon- 
struction. 

We have verified our calculations against results ob- 
tained with CMBFAST version 2.4.1 for models sup- 
ported by the latter (open and flat). We have also com- 
pared our results with a pre-release version of CMBFAST 
which Seljak & Zaldarriaga are developing to support 
closed models, with good agreement well into the damping 
tail. 

3. RESULTS 

In Figure |^ we plot the intensity and polarization power 
spectra in ACDM models assuming no reionization. One 
model is closed (i^tot = 1-2, f^A = 0.8), while the other is 
flat (r^tot — 1, f^A — 0.6). In both cases we take the mat- 
ter fraction = 0.4, baryon fraction fib = 0.045, and 
Hubble's constant Hq = 65kms~^Mpc~"'^. For the scalar 
modes, we assume adiabatic initial conditions with a scale- 
invariant primordial power spectrum 7^,s(t^) = constant. 



Our c onventions generalize those of Lyth fc Woszczynq 



(1995 ), so that the gradient of the 3-Ricci scalar on co- 
moving hypersurfaces DaR^^^ , receives power 



u>3 



1/^-4 



(5) 



from super-Hubble modes. The tildes denote that the 
quantity on the left is evaluated in the (energy) frame 
where the momentum density qa vanishes. For the ten- 
sor modes, we assume a scale-invariant spectrum 'Pt{v) = 
constant. Our conventions are such that the power in the 
electr ic part of the Weyl tensor Eab (e.g. Ellis fc van eIs^ 
1999 ) from super-Hubble modes is 



{EabE' 



4) 



(6) 

Note that in closed models, the infra-red divergence in 
the tensor power spectrum seen in open models with 
V'Y(y) = constant is avoided because of the cutoff a.tv = 'i. 

The models in Figure |] have equal physical densities, 
r^m-ffo 3''^'^ f^b-ffoi so the sound horizon at last scatter- 
ing and the early-time dynamics are approximately equal 
in the two models. On small angular scales, where the 



polarization and scalar anisotropies are projections of ef- 
fects at last scattering, restricting parameter changes to 
r^A and the curvature leads to approximate scaling of the 
CMB power spectra: Ci Cat where a is the ratio of the 
angular-diameter distances to last scattering in the orig- 
inal and final models. Approxim ations of this sort were 



used by Efstathiou et al. (1999 ) to include closed mod 
els in their recent joint analysis of CMB and supernovae 
data. On large angular scales (small ^), the approximate 
scaling is broken for scalar anisotr opies by the (late-time ) 
integrated Sachs- Wolfe effect (e.g. Hu fc Sugiyama 1995 ), 
which is responsible for the enh ancement at low t seen in 
the closed model in Figure |l| ( [Abbott fc Schaefer 1986 ; 
White fc Scott 1996). The effects of curvature terms in 



the primordial power spectrum arc a lso potentially observ- 
able on large scales ( Hu et al. 1998| ). The presence of the 
curvature scale in the primordial power spectrum, com- 
bined with the cutoff at v ~ 3, suppresses the intensity 
quadrupole for tensor modes in the closed model relative 
to the flat model in Figure [l](c.f. the case of open models, 
e.g. [Hu fc White 199 < 




Fig. 1. — Scalar (left) and tensor (right) intensity and polarization 
power spectra in a closed CDM model (Ctot = 1-2, = 0.8; thin 
lines), and a fiat model (Htot = 1, f^A = 0.6; thick lines). In both 
cases Hm = 0.4, = 0.045, and _ffo = 65kms~'^Mpc~^. The upper 
solid lines are the intensity, the lower ones the electric component of 
the polarization C^^ , and the dashed lines the magnetic component 
C^^ . (Note scalar modes produce no magnetic polarization.) The 
scalar and tensor intensities are normalized to unity at £ = 10. 



4. CONCLUSION 

We have presented the first calculation of the CMB 
power spectra, including the effects of polarization, in 
closed FRW models. We have implemented the efficient 
line of sight algorithm for general geometries, using co- 
variantly defined, gauge-invariant variables, thus allow- 
ing accurate and rapid modelling over the full volume of 
parameter space of FRW models. Our Fortran 90 code, 
ba sed on CMBFAST version 2.4.1, is publically available 
at ^rttp://www. mrao.cam.ac.uk/~amll005/cmb . 
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